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ABSTRACT 

We examine the dynamics structure of the rich cluster A1689, combining VLT/VIMOS spectroscopy 
with Subaru/Suprime-Cam imaging. The radial velocity distribution of ^ 500 cluster members is 
bounded by a pair of clearly defined velocity caustics, with a maximum amplitude of ~ |4000| km/s at 
~ 300 h~^ kpc, beyond which the amplitude steadily declines, approaching zero velocity at a limiting 
radius of ~ 2 h~^ Mpc. We derive the 3D velocity anisotropy and galaxy number density profiles using 
a model-independent method to solve the Jeans equation, simultaneously incorporating the observed 
velocity dispersion profile, the galaxy counts from deep Subaru imaging, and our previously derived 
cluster mass profile from a joint lensing and X-ray analysis. The velocity anisotropy is found to be 
predominantly radial at large radius, becoming increasingly tangential towards the center, in accord 
with expectations. We also analyze the galaxy data independently of our previous analysis using two 
different methods: The first is based on a solution of the Jeans equation assuming an NEW form 
for the mass distribution, whereas in the second method the caustic amplitude is used to determine 
the escape velocity. The cluster virial mass derived by both of these dynamical methods is in good 
agreement with results from our earlier lensing and X-ray analysis. We also confirm the high NEW 
concentration parameter, with results from both methods combined to yield Cvir > 13 (Ic). The 
inferred virial radius is consistent with the limiting radius where the caustics approach zero velocity 
and where the counts of cluster members drop off, suggesting that infall onto A1689 is currently not 
significant. 

Subject headings: clusters: A1689 - clusters: dynamical analysis - clusters: galaxies, DM 
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1. INTRODUCTION 

Clusters of galaxies display a wide range of distinct 
observational phenomena, providing detailed physical 
information of central importance to cosmology. Re- 
cent observations have resulted in an impressive consis- 
tency with the predictions of the standard ACDM model 
(Tegmark et al. 2004; Spergel et al. 2007 and references 
therein). In this model, a cosmological constant domi- 
nates the cosmic energy budget today, but galaxies and 
other structures were assembled earlier, primarily out of 
cold dark matter (CDM). While the model successfully 
matches observations of the primary anisotropy of the 
cosmic microwave background and the large-scale struc- 
ture in galaxy surveys, it is also important to test its 
validity on smaller scales. The abundance and structure 
of non-linear objects are sensitive probes of the prop- 
erties of dark matter (DM) and of the primordial den- 
sity fluctuation field. Gas cooling and various feedback 
mechanisms complicate the interpretation of the DM dis- 
tribution in galaxies (e.g.. Bower et al. 2006). However, 
for clusters most of the gas is observed to be too hot 
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and rarefied to cool efficiently and is thus expected to 
be in hydrodynamical equilibrium at the virial tempera- 
ture, and therefore to have only a small influence on the 
non-baryonic mass distribution (outside the inner core). 
Hence it is reasonable to suppose that constraints on the 
mass profile of a cluster probe the dominant DM. 

Interaction of a cluster with another cluster or group of 
galaxies is commonly indicated by the presence of shock 
fronts, seen in high quality X-ray observations. In the 
case of the "bullet cluster" , a cone-shaped shock front is 
visible indicating that two clusters have passed through 
each other, with an obvious merger of the intracluster 
(IC) gas of the two clusters, but the galaxies and the 
lensing mass distribution are largely intact. This di- 
rectly implies that the DM is coUisionless like the galax- 
ies (Markevitch et al. 2002; Clowe et al. 2004; Bradac 
et al. 2006), supporting the simplest possibility that DM 
interacts only via gravity. 

For relaxed clusters, knowledge of the mass profile can 
provide crucial information regarding the nature of DM 
and the thermal history of IC gas. Precise measurements 
have been made by dynamical analysis of galaxy veloc- 
ities, via hydrostatic analysis of X-ray observations, or 
directly via lensing. Detailed measurements of the X-ray 
spectrum and intensity profile yield the distribution of 
total mass through solving the equation of hydrostatic 
equilibrium. In practice the temperature measurements 
are often made uncertain by significant complications 
such as the likely multiphase nature of gas in equilib- 
rium. We have recently developed a model-independent 
joint lensing/X-ray analysis (Lemze et al. 2008, hereafter 
L08) to examine the consistency of the X-ray tempera- 
ture and surface brightness profiles with the lensing data. 
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finding that the cluster's spectraUy measured tempera- 
ture profile is systematically ~ 30 — 40% lower than de- 
duced from solving the equation of hydrostatic equilib- 
rium with the precisely measured lensing mass profile. 
This may reflect in part the ambiguity in deriving 3D 
temperatures from projected X-ray data, since a given 
line of sight will in general intersect a range of gas tem- 
peratures if the gas is not strictly isothermal (Mazzotta 
et al. 2004; Vikhlinin 2006). It is also conceivable that 
the gas is not strictly single phase but may contain small 
scale structure, including relatively dense cooler clouds 
as found in detailed simulations (Kawahara et al. 2007), 
which may lead to a significant downward bias in ob- 
served temperature estimates. 

The velocity dispersion profile of clusters has long been 
used to estimate the galaxy dynamics and the cluster 
mass profile via the Jeans equation, when sufficient red- 
shift information is available. This technique was orig- 
inally applied to several well studied nearby clusters 
(Fuchs & Materne 1982; Sharpies et al. 1988), and more 
recently to the results of dedicated surveys (Carlberg et 
al. 1997; Katgert et al. 2004; Biviano & Katgert 2004; 
Hwang & Lee 2008). In these studies the expected veloc- 
ity anisotropy complicates the interpretation, as orbits 
are not expected to be isotropic but to become predom- 
inantly radial towards the cluster's virial radius. More 
recently, with the availability of larger samples of redshift 
measurements, it has been recognized that relaxed clus- 
ters should have sufficiently well defined velocity "caus- 
tics" , providing an independent means of deriving clus- 
ter mass profiles; this has been applied to several clusters 
and compared with lensing-based masses (Geller, Diafe- 
rio, & Kurtz 1999; Rines et al. 2003; Diaferio, Geller,& 
Rines 2005). Joint studies of galaxy dynamics and lens- 
ing in the central cluster region are also providing de- 
tailed mappings of the central mass distributions of mas- 
sive clusters, where the velocity dispersion profile of the 
central galaxy can be compared with the analysis of mul- 
tiply imaged sources (Sand et al. 2008). 

Lensing work is now able to measure cluster mass pro- 
files with sufficient precision to usefully test the distinc- 
tive prediction of a relatively shallow mass profile for 
halos dominated by CDM (Navarro, Frenk, & White 
1997, hereafter NEW; Hennawi et al. 2007; Duffy et 
al. 2008). Combined weak and strong lensing measure- 
ments have shown that the continuously steepening form 
of the NFW profile is a reasonable description for the 
mass profiles of three carefully studied clusters (Kneib 
et al. 2003; Gavazzi et al. 2003; Broadhurst et al. 2005a, 
hereafter B05a), although with surprisingly high values 
derived for the profile's characteristic concentration pa- 
rameter in each case (B05a). More recently, the concor- 
dance ACDM model has been examined with increasing 
(though still relatively small) samples of clusters, indi- 
cating that the concentrations derived are significantly 
higher than predicted over a wide range of cluster mass, 
after a statistical correction for lensing-induced biases. 
This is seen both in terms of the size of the Einstein ra- 
dius (Broadhurst & Barkana 2008) and the weak lensing 
profiles of several well known clusters (Broadhurst et al. 
2008). This question has also been examined recently 
by stacking the lensing signal of cluster samples identi- 
fied in the SDSS survey. A wide dispersion in derived 
concentrations is found, with an average value in better 



agreement with the predicted relation (Mandelbaum et 
al. 2008). Here it is crucial that reliance on photometric 
redshifts does not result in the misclassification of cluster 
members as background galaxies, which would artificially 
decrease the central lensing signal thereby lowering mea- 
sured concentrations. 

Here we concentrate on the internal dynamics of 
A1689, one of the best-studied clusters, using a relatively 
large redshift survey from the VLT/VIMOS wide- field in- 
strument. This cluster is attractive as it is very massive 
and appears to be relaxed, with a highly symmetric X- 
ray morphology (L08; Riemer-Sorensen et al. 2008, see 
their fig. 1), the centroid of which coincides both with 
the cD galaxy and the center of mass derived from lens- 
ing (Xue & Wu 2002; L08). Only a low level of sub- 
structure is visible in detailed lensing and X-ray maps 
(B05a; Broadhurst et al. 2005b; L08; Umetsu & Broad- 
hurst 2008), though Andersson & Madejski (2004) claim 
an indication for a small deviation from a relaxed state. 
Earlier dynamical analyses were made by den Hartog & 
Katgert (1996), Czoske (2004), and Lokas et al. (2006). 
In our earlier work we found that the cluster has a rela- 
tively high concentration parameter, which makes this 
a very interesting target to explore with independent 
dynamical means. We also make use of the extensive 
VLT/VIMOS data to explore galaxy dynamical proper- 
ties in unprecedented detail, including in particular the 
velocity anisotropy profile, which has been difficult to de- 
termine from previous observations (Benatov et al. 2006). 

This paper is organized as follows. In §[2] we describe 
the data, from which wc derive the radial profile of the 
projected surface density (ii l2.1|) . and the projected veloc- 
ity distribution and radial profile of the velocity caustics 
(§ 12. 2p . In § [3] we detail our method of using the Jeans 
equation together with previous knowledge of the clus- 
ter mass profile to determine the three-dimensional radial 
profiles of galaxy velocity dispersion and anisotropy, with 
results given in § 2) In § O we independently determine 
the cluster mass profile from the data on galaxy dynam- 
ics, using the caustics as a measure of the escape velocity 
from the cluster 15. 1|) . or fitting an NEW profile to the 
dynamical data fi) l5.2p . We compare in §[6] various defini- 
tions of the cluster boundary, based on the radial profiles 
of the cluster mass, galaxy density, or caustics, and end 
with a discussion in § [T] 

2. DATA ANALYSIS 

We use several different types of data in this pa- 
per: galaxy positions and velocities, gravitational lens- 
ing, both strong and weak, and X-ray emission and spec- 
troscopy. Previously we have made joint use of the same 
lensing and X-ray data to derive a cluster mass profile by 
a model-independent joint analysis as described in L08. 
Here we extend this work and analyze independent dy- 
namical information on the galaxy velocities and also the 
number density profile of cluster members. In this sec- 
tion we explain in detail our analysis of the latter two 
data sets. 

2.1. Galaxy surface number density 

Establishing the form of the projected profile of clus- 
ter member galaxies is essential for fully employing the 
Jeans equation (Binney & Tremaine 1987; see §[3]), which 
allows for a spatial distribution of galaxies that need not 
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follow the dominant DM. Measurement of the projected 
galaxy distribution requires subtraction of the back- 
ground and foreground field galaxy populations. This 
must be achieved with accurate multi-color photometry 
as spectroscopy is usually limited to small samples and 
does not extend faint enough to include the majority of 
cluster members. 

We used Subaru photometry in the V and I bands. The 
data are complete to a depth of Iab = 26.5, reaching 
eight magnitudes below L* of the luminosity function. 
We selected a color region based on a color-magnitude 
relation analysis, which comprises 17474 galaxies. This 
color region ranges from the red side of the E/SO cluster 
sequence to a blue boundary chosen so that the sample 
extends sufficiently blueward to include most (~ 85%) 
of the cluster members as described in Medezinski et al. 
(2007, hereafter M07, fig. 1). M07 established that this 
color region contains the majority of cluster members in 
addition to some background galaxies, by examining the 
weak lensing signal. In this color region, which includes 
the E/SO sequence and bluer objects (including some 
background galaxies) the lensing signal was found to be 
significantly lower than the true background signal as 
measured for red background galaxies. This is because 
(unlensed) cluster members dilute the weak lensing signal 
of the background, allowing us to identify the region of 
color space occupied by cluster members. 

We derived the projected galaxy distribution of clus- 
ter members from the above color-selected galaxy sample 
based on fitting the radial distribution. Galaxies were 
first radially binned into annuli and the surface number 
density in each annulus was determined; note that at 
large radii only part of each annulus was covered by the 
detector. Assuming Poissonian (V^) errors, the galaxy 
surface number density was modelled as that due to the 
cluster plus a background galaxy surface number den- 
sity, with the latter assumed uniform. Thus, we effec- 
tively neglected fluctuations in the background galaxy 
surface number density on scales smaller the cluster size 
(about 22' at z = 0.183). There could be such fluctua- 
tions at some level due to correlated structures (such as 
filaments) along the line of sight to the cluster. Galax- 
ies in these structures could mistakenly be included as 
cluster members, a possibility that we assess below. 

Figure [1] shows that the radial profile of the above 
color-selected galaxy sample is well fitted with a uniform 
background plus a general cored profile to represent the 
cluster galaxy surface density profile: 

^0 

Stot = Sgai + Cbg = j-^ ^ (r/rc)2]P ^''^ ' 

where Etot is the total galaxy surface number density, Eq, 
Tc, and p are the three parameters of the cored profile, 
and Cbg is the background density. The resulting fit is 
good, x^/dof = 19.2/16. Using a larger number of bins 
than A*" = 20 improves the reduced , e.g., = 50 
and A^ = 100 give xV^of = 36/46 and xVdof = 96/96, 
respectively. The level of the background is stable and 
does not depend on the number of bins, e.g., for A^ — 20, 
50, and 100 we find C = 1071 ± 40, 1064 ± 42, and 
1063 ±41, respectively. We settled on just 20 bins, since 
this made it easier to combine the surface density data 
set with the projected velocity dispersion data set, which 
is inferior in terms of signal to noise per radial bin (see 



below). With 20 bins, we found best-fit values of Eq = 
1200 ± 140 h^ Mpc-2, Tc = 450 ± 150 h^^ kpc, and 
p = 0.74 ±0.30. 
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Fig. 1. — Galaxy surface number density. We show the to- 
tal (blue circles), the cluster contribution (red triangles), and the 
background value (solid horizontal line) and its uncertainty (dashed 
lines). All uncertainties are 1-cr. We also show the core profile that 
best fits the cluster galaxy surface number density (dash-dotted 
curve), and the corresponding fit to the total galaxy surface num- 
ber density (dashed curve). 



2.2. Projected Velocity Dispersion 

To measure the mass profile of galaxy clusters using 
galaxy motions it is necessary to obtain precise velocity 
measurements for a statistically large sample of galaxies. 
The data used here are part of an extensive multi-object 
spectroscopy survey carried out with the VIMOS spec- 
trograph on the VLT (Czoske 2004); for observational 
details, see Czoske (2004). This dataset constitutes 1469 
objects with reliable spectroscopic redshifts, a major ad- 
vance over previous surveys of A1689. Note though that 
we did not try to find the galaxy surface density pro- 
file from the projected velocity data set, since the data 
set used in the previous section contains a much larger 
number of galaxies. 

We analyzed the spectroscopic sample by first defining 
cluster membership using the velocity "caustics" , which 
are clearly visible for this cluster in the form of a bound- 
ary which varies with radius, peaking at around ±4000 
km/s at ~ 300 h~^ kpc and declining steadily at larger 
radius (see figured]). The caustics are related to the es- 
cape velocity from the cluster and thus provide a tangible 
physical basis by which we can separate cluster mem- 
bers from foreground and background galaxies. Defining 
membership is especially important for massive clusters 
like A1689, since they have a relatively wide internal ve- 
locity spread and they extend to large radii; this increases 
the chance for interlopers, which can have a large effect 
on the derived projected velocity dispersion, especially at 
large radii where the number density of cluster members 
is low. Wojtak & Lokas (2007) showed that without the 
proper removal of interlopers the inferred parameters of 
the mass distribution in the cluster are strongly biased 
towards higher mass and lower concentration. See also 
Wojtak et al. (2007) for interesting comparison between 
different methods of interlopers removal, which unfortu- 
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nately do not include the D99 caustics approach used in 
this paper. 

To define the caustics we employed the technique pio- 
neered by Diaferio (1999, hereafter D99) based on a mul- 
tidimensional adaptive kernel method (Silverman 1986; 
Pisani 1993; Pisani 1996). A summary of the applica- 
tion of this technique can be found in Diaferio, Geller, 
& Rines (2005); it has previously been applied to sev- 
eral clusters (Reisenegger et al. 2000; Biviano & Girardi 
2003; Diaferio, Geller, & Rines 2005; Rines & Diaferio 
2006). Briefly, a redshift-space diagram is constructed, 
i.e., the line of sight velocity v is plotted versus the pro- 
jected distance R from the cluster center, and for a well 
defined cluster these points should be distributed in a 
characteristic "trumpet" shape, the boundaries of which 
are termed caustics (Kaiser 1987; Regos & Geller 1989). 
The D99 procedure locates the caustics and determines 
the radial dependence of their amplitude (in units of ve- 
locity), which is related to the escape velocity and thus 
depends on the mass profile. Galaxies that are inside the 
caustics can then be considered to be cluster members. 
Of course, some of these galaxies might be interlopers, 
but their number is typically a few percent and has little 
effect on dynamical analyses (A. Diaferio, private com- 
munication) . 

To apply this procedure we determined the threshold 
K that defines the caustic location through fq {R, v) ^ k 
(D99). Here fq{R,v) is the galaxy density distribution 
in the redshift-space diagram, smoothed with an adap- 
tive kernel, and v is the peculiar velocity (measured with 
respect to the cluster's mean redshift). The parameter 
q sets the scaling between the quantities R and v within 
the smoothing procedure. We used q — 27, close to the 
value usually used, 25 (D99; Rines et al. 2003). We note 
that different values of q in the 10 — 50 range have little 
effect on the results (D99). 

The parameter k was chosen by minimizing the quan- 
tity S{k, (R)) = |(i;L).,(fl)-4(i;2)|2, where {vI,),,j^r) = 

J^"^ A^{R)(l){R)dR/ J^'^^ (l}{R)dR is the mean value of 
the square of the caustic amplitude A{R) within (i?) 
(the mean projected radius of the cluster members), 
(/)(i?) = / fg{R,v)dv, and (u^) is the one-dimensional 
velocity dispersion of the cluster members. The uncer- 
tainty of A is proportional to the inverse of the galaxy 
number density within the caustics, and was estimated as 
SA(R)/A{R) — k/ max{/g(i?, w)}, where the maximum 
is found along the -y-axis at each R (D99). 

Figure [2] shows the velocity-space diagram and radial 
profile of the derived velocity caustics. There are a total 
of 476^43 galaxies identified as cluster members lying be- 
tween these caustics. The velocity-space diagram clearly 
looks qualitatively different inside and outside the caus- 
tics, indicating that the caustic-finding procedure has 
identified a physically-meaningful boundary. The den- 
sity of the galaxies that lie between the caustics shows 
a smooth decline with distance from the center, indicat- 
ing that the cluster has a simple monolithic dynamical 
structure. The radius at which the caustics meet at zero 
velocity represents an effective maximum radius for the 
cluster, which physically includes cluster galaxies both 
within the virial radius and others that may be currently 
infalling onto the cluster (see the discussion in § [7]) . Note 
that the two caustics are by definition symmetric about 



u = in this method, which figure [2] clearly shows is a 
valid assumption for this cluster . 




Fig. 2. — Velocity-space diagram of A1689. The caustics are 
shown along with l-c error bars (shown only on the outer side of 
the caustics for clarity). 



Having defined cluster membership based on the loca- 
tion of the caustics, we can estimate the projected veloc- 
ity dispersion profile restricted to cluster members. We 
divided the cluster members into velocity bins of size 600 
km/sec, sufficiently large compared to the redshift mea- 
surement uncertainty of ~ 200 km/sec. Different bin- 
ning gave very similar results. We fitted a Gaussian to 
the velocity histogram to calculate the overall projected 
velocity dispersion (figure [3]); it provided a reasonable fit 
with a reduced (hereafter Xr) of 1.7 per degree of free- 
dom. If we ignore the caustics and simply bin all galaxies 
within the maximal velocity range \v\ < 4000 km/sec, the 
fit to a Gaussian is significantly worse, Xr — 2-5, with a 
more asymmetric distribution and a somewhat larger ve- 
locity dispersion. It is important to mention that we did 
not use the results of the Gaussian fit in our analysis. 
We show this fit to argue in favor of the real physical 
meaning of the caustics as we have measured them, and 
demonstrate that without excluding unrelated interlop- 
ers with large velocities, dynamically derived quantities 
of the cluster may be significantly skewed, particularly 
at large projected radii. Note that we fitted a Gaus- 
sian for simplicity, although the real distribution may be 
slightly non-Gaussian (see, e.g., Kazantzidis, Magorrian 
& Moore 2004; Sanchis, Lokas & Mamon 2004; Diemand, 
Moore & Stadel 2004; Wojtak et al. 2005). 

We next estimated the radial profile of the projected 
velocity dispersion from our sample of cluster members 
identified within the caustic boundaries. We divided the 
data into 10 radial bins, so that in each bin there were 
at least 47 galaxies. The radius assigned to each bin was 
the mean of the measured radii of galaxies lying within 
the bin. The projected velocity dispersion was taken 
to be the standard deviation of the galaxy velocities in 
the bin. The error in the projected velocity dispersion 
Aap{R) was calculated as {Aap{R))'^ = (A„iofTp(-R))^ -I- 

{Ass,ap{R)f + {AccCrp{R)f , where A^^apiR), Asi,ap{R), 
and AccO'p(-R) are the measurement, sample, and cluster 
membership uncertainties, respectively. The measure- 
ment error is A„ieO-p(_R) — l/ap{R) ■ dvi/{N{R) — 1), 



5 



120 




V [km/sec] 



Fig. 3. — Observed histogram of galaxy velocities and Gaussian 
fits. Top panel: all galaxies; the fit yielded Xr — 2.5 with derived 
parameters Zclustcr = 0.185 ± 0.0003 and ap = (1.75 ± 0.08) X 10^ 
km/sec. Bottom panel: only cluster members, whose identification 
is based on the caustics; the fit yielded Xr — ^-"^ with Cp = (1.40 ± 
0.06) X 10^ km/sec. 

where N{R) is the number of galaxies in the bin, and 
dvi = 200 km/sec is the measurement error of the pro- 
jected velocity of the ith galaxy. The 1-cr sampling error 

is A,,ap±{R)/ap{R) = [(1-2/(9z/)T\/27(9^)"'/'-1], 
where this formula is valid for i/ > 30 with — N{R) — 1 
the number of degrees of freedom (Danese, De Zotti, & 
Tullio 1980). AccCrp{R) is the error in the projected ve- 
locity dispersion from errors in estimating which galax- 
ies are cluster members and which are interlopers, aris- 
ing from the uncertainty in the amplitude of the caustic, 
SA{R). We estimated Acc(Jp{R) by calculating the differ- 
ence of the projected velocity dispersion when the clus- 
ter members are defined to lie within the caustic defined 
by the amplitude A + 6 A instead of A: AccTp ±{R) — 
\apiA±SA)-ap{A)\. ^ 

There are two previous published measurements of 
the projected velocity dispersion for A1689. Lokas et 
al. (2006) took all galaxies in the NASA Extragalactic 
Database (NED) with redshifts z = 0.1832 ± 0.05 and 
located at projected distances below 2 Mpc from the 
cluster center. The redshift data for the galaxies came 
from different surveys, mainly by Teague, Carter, & Gray 
(1990), Balogh et al. (2002), and Due et al. (2002). While 
the sample of Teague, Carter, & Gray (1990) came from 
a standard magnitude-limited survey, the selection cri- 



teria of those of Balogh et al. (2002), and Due et al. 
(2002) were aimed at star-forming galaxies, which may 
bias the sample towards outer regions with more sub- 
structure. This full sample comprised 192 galaxies, which 
they reduced to 130 after defining a constant velocity 
cutoff of 3000 km/s, estimated visually, which is sub- 
stantially smaller than the 4000 km/s maximum ampli- 
tude of our caustics. On this basis they claimed signifi- 
cant substructure in velocity along the line of sight, but 
there is no evidence for such substructure in our velocity 
data (Czoske 2004), nor is significant substructure seen 
in the X-ray emission maps or the lensing-based mass 
map (L08; Umetsu & Broadhurst 2008). In an earlier 
work den Hartog & Katgert (1996) determined the ra- 
dial velocity profiles of 72 clusters, including A1689. For 
A1689 they used the data from Teague, Carter, & Gray 
(1990) comprising only 63 cluster members. The pro- 
jected velocity profiles derived in these two papers are 
compared to the one we derived in figure IH 

3. METHODOLOGY 

In this section we present our procedure for explor- 
ing the structure of the cluster using the above veloc- 
ity data and the projected profile of the galaxy distri- 
bution derived above. In addition, we make use of the 
mass distribution that we derived in earlier work from a 
combined X-ray and lensing analysis (L08). In L08 we 
combined lensing and X-ray measurements to determine 
model independent profiles of the total mass and the gas 
mass. Here we add the data on the galaxy surface num- 
ber density and the projected velocity dispersion, and 
relate them via the Jeans equation. From this we obtain 
the 3D galaxy number density profile and the velocity 
anisotropy profile. Our work is the first model-free de- 
termination of the galaxy velocity anisotropy in a cluster. 

In general, we determine the best-fit values of our free 
parameters by fitting the profiles of galaxy surface num- 
ber density and projected velocity dispersion simultane- 
ously. We first proceed with simple analytical forms for 
the relevant profiles, in which we fit for the free param- 
eters of these expressions. We then proceed to a more 
flexible model-independent approach, similar to that de- 
veloped in L08. In our model-independent approach, the 
free parameters are simply the values of the 3D profile 
of the galaxy number density and velocity anisotropy at 
several fixed (equally spaced) radii. The radial ranges of 
these free parameters are set according to the total span 
of each data set. The values of each quantity (i.e., the 
galaxy number density or velocity anisotropy) across the 
cluster are given by linearly interpolating from its val- 
ues at the fixed radii. Thus, the method only assumes 
smoothness. 

The galaxy number density and velocity dispersion are 
related by the Jeans equation for a steady-state spheri- 
cally symmetric system, 

dr r 

(2) 

where ngai(r) is the galaxy number density, and P{r) is 
the velocity anisotropy parameter: 

/3(r)^l-474' (3) 
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where crr{r) and crt(r) are the radial and the tangen- 
tial components of the velocity dispersion (Binney & 
Tremaine 1987). The measured quantities used here are 
the galaxy surface number density, Sgai, 



Sgal(i?) = 2 



ngi,i{r)rdr 



i?2 



(4) 



and the observed projected velocity dispersion, 



Sgal(i?) 



»^gaiWo-^(r) l-/3(r) 



rdr 



(5) 

Given a total cluster mass profile M{r), for any as- 
sumed profiles of rigai and (3 we use the Jeans equation 
(eq. ^ to derive the profile of the radial velocity dis- 
persion ijr- We then use eqs. Q and ([5]) to compare 
the model to the two observed data sets, and find the 
best-fit free parameters. This procedure was first ap- 
plied using simple analytic forms for the profiles of 3D 
galaxy number density ngai and velocity anisotropy pa- 
rameter /3. The model-independent mass profile found 
in our earlier work (LOS) was used; we also compared in 
some cases the effect of using its approximation as an 
NFW profile, with Cvir = 12.21?-^ and w,, = 2.14t°J^ 
h~^ Mpc for the concentration parameter and virial ra- 
dius, respectively (LOS). For the galaxy number density 
profile we assume the (3 model 



r^gal(r) = ?T.gal(0) [l -I- {r/vsY 



(6) 



where here we use p instead of the standard (3 (which 
we use for the velocity anisotropy parameter); e.g., the 
King profile (King 1962) corresponds io p = 1. N-body 
simulations for a variety of cosmologies show that the 
velocity anisotropy has a nearly universal radial profile 
(Cole & Lacey 1996). We follow Carlberg et al. (1997) 
and model it as: 



(C+1) 



(r/rvc)' 



(r/r 



1 



C 



(7) 



where C and r^c (a core radius of the velocity profile) 
are free parameters. The free parameter C is needed 
to ensure the freedom for f3 to take on any value in its 
allowed range, — oo < /3 < 1, including circular orbits at 
small radii. A similar profile is advocated by Girardi et 

al. (199S) for "case A" clusters, (3 = ^.-Z^^i , except that 
we do not constrain the velocities to be isotropic in the 
center. 

When performing model-independent fits, we must ex- 
trapolate each quantity beyond the last data point. We 
extrapolated the anisotropy parameter /3 to be constant, 
i.e., equal to its value at the last projected velocity dis- 
persion data point. In this way we did not pre-constrain 
it to equal unity at large radii. We extrapolated the 



galaxy number density oc 



consistent with a King 



profile and also with the asymptotic behavior of our best- 
fit f3 model (which is cx 7.-3 i8±o.42. table [1] below) . 

4. THE DYNAMICAL PROPERTIES OF THE 
CLUSTER GALAXIES 

In this section we present the results for the dynami- 
cal properties of the galaxies, namely the galaxy number 



density, velocity anisotropy, and radial velocity disper- 
sion, all derived from our simultaneous analysis of the 
galaxy surface number density and velocity dispersion 
data sets. As stated in the previous section, we em- 
ployed the Jeans equation and used the mass profile we 
had previously derived from lensing and X-ray data (with 
our model- independent method described in LOS). We 
present here the results both from fitting simple analyt- 
ical forms to the dynamical data and from our model- 
independent approach. 

Using the assumed profiles, the simultaneous fit to 
the two data sets was very good, as expressed by the 
low Xr = 15.6/(30 — 5), where there are 30 total data 
points and 5 parameters (from eqs. [B] and [7]). The 
contribution of each data set within the full simulta- 
neous fit to both, was Xr('^uniber) = 12/(20 — 3) and 

(velocities) = 3.6/(10 — 2). The fits are shown in 
figure m The best-fit parameters of the analytical ex- 
pressions for the galaxy number density and velocity 
anisotropy profiles are given in table [1] 
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den Hartog & Katgert data 
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Radius [fi^^ l^pc] 

Fig. 4. — Data and best fits based on analytical profiles. Top 
panel: The observed profile of galaxy surface number density (blue 
circles) is compared with our best-fit model profile (red curve). 
Bottom panel: The observed profile of projected velocity disper- 
sion (blue circles) is compared with our best-fit model profile (red 
curve). We also show the measurements of Lokas et al. (2006) 
(black triangles) and den Hartog & Katgert (1996) (green squares). 
In both panels, l-cr measurement uncertainties are indicated. 



The derived 3D profile of galaxy number density is 
shown in figure [5] The analytical cored profile is con- 
sistent with our best-fit model-independent profile. We 
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TABLE 1 

Best-fit parameters for the galaxy number density and 

VELOCITY ANISOTROPY PROFILES 



Parameter 


Value 


ngal(O) [h« Mpc-»] 
Ts [h~^ kpc] 
P 

rvc [h~^ kpc] 
C 


1380 ± 320 
455 ± 110 
1.06 ±0.14 
395 ± 280 
1.8 ±2.6 



also tried a cuspy profile as an analytical form for the 
galaxy number density, instead of a cored profile. Specif- 
ically, we adopted a general form with a l/r cusp, 
"■gai(r) = ngai(0)/[(r/rs)(l -|- {r/rs))P], which includes 
NFW [p = 2) and Hernquist (1990) [p = 3) profiles as 
special cases. Although this profile gave an acceptable 
fit, Xr = 26.3/(30 — 5), the fit was worse than that of the 
cored profile. The best-fit model gave ngai{0) = 18.8 and 
an unrealistic value for Ts, Ts = 15 h~^ Mpc; with this 
value of Ts and the deduced p, p — 18.1, the p-dependent 
term is negligible (implying that only the ratio Ugai (0) /r^ 
is relevant). 




500 1000 1500 2000 

Radius [h^^ kpc] 



Fig. 6. — Reconstructed 3D velocity anisotropy profile. We show 
the best-fit result and l-a range from fitting the simple profile of 
eq. Q (solid and dashed curves), where the best-fit model parame- 
ters are listed in Table [T] Also shown are the values and l-cr errors 
as derived from the model-independent approach, where we used 5 
free parameters for the values of the galaxy number density and 5 
(red circles) free parameters for the values of /3 at 5 equally spaced 
radii. The value of the velocity anisotropy was extrapolated as 
constant beyond the last data point. 




10^ 10^ 10* 

Radius [h^^ l^pc] 



Fig. 5. — Reconstructed 3D galaxy number density profile. We 
show the best-fit result and l-a range from fitting a /3 model (solid 
and dashed curves), where the best-fit model parameters are listed 
in Table[T] Also shown are the values and l-cr errors as derived from 
the model-independent approach (red circles), where the free pa- 
rameters were the values of the galaxy number density at 5 equally 
(linearly) spaced radii. 

The velocity anisotropy profile /3 is shown in figure [HI 
Again, the result from fitting a particular simple ana- 
lytical form for the profile agrees with the result of the 
model-independent approach which does not assume any 
particular shape for the profile. While the uncertainties 
are large, the figure indicates a clear tendency of an in- 
creasing P with radius. The ct^ profile, as obtained from 
the Jeans equation (eq. H]) , is shown in figure [T] 

The analytical expression we used for /3, eq. ([7]), con- 
strains f3 to equal unity at large radii, corresponding to 
purely radial orbits. D99 analyzed N-body simulations 
and derived values somewhat smaller than 1 at their lim- 
iting radius, 6r2oo- To allow for deviation from unity we 

also fit the expression (C + a)j^^j{rjT^ " C, where a is 
a free parameter which governs the asymptotic behavior 



io' 




Radius [h ^ l^pc] 

Fig. 7. — Reconstructed profile of the radial velocity dispersion 
fTr- We show the best-fit values and l-cr uncertainties from fitting 
with the simple analytical expressions in eqs. JHJ and jTj (solid and 
dashed curves). 

of (3. The value of a = 1 gave the best fit, also consis- 
tent with the value obtained by the model-independent 
method at the largest radial point. 

5. DERIVING THE MASS PROFILE 

In § [3] we used the mass profile as derived from X-ray 
and lensing data together with the new data on the pro- 
jected velocity dispersion and the galaxy surface number 
density to derive the 3D profiles of the velocity anisotropy 
and galaxy number density. Alternatively, there are at 
least two different ways to derive the total mass profile 
of the cluster directly from the data sets on galaxy dy- 
namics. Perhaps the simplest approach is to use only the 
velocity caustics derived above, the amplitude of which 
is related to the escape velocity, which is a tracer of the 
cluster mass profile. The second approach is to use both 
data sets (of the projected velocity dispersion and galaxy 
surface number density) together with the Jeans equation 
to fit an NFW mass profile. Thus, we have two indepen- 
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dent methods for estimating the cluster mass profile. 

5.1. Mass profile from velocity caustics 

As mentioned already, D99 has shown that the 3D mass 
profile can be fairly well estimated directly from the am- 
plitude of the velocity caustics. In particular, the total 
mass within radius r is estimated as 



M{< r) = ^ J\\R)dR 



(8) 



This equation essentially estimates the mass profile based 
on the local escape velocity, but note that the prefactor of 
0.5 is not the result of an exact derivation. This prefactor 
comes from the more general suggestion by Diaferio & 
Geller (1997; hereafter DG97) for the relation between 
A and M: M(< r) = ^ A'^{R)FpiR)dR in terms of 



p(r)r2 3 - 2/3(r) 
l-/3(r) 



(9) 



where p(r) is the cluster mass density and (f)(r) = 
_GA^<r) _ p{x)xdx is the gravitational poten- 

tial generated by the cluster. D99 noted that the func- 
tion Fpir) is slowly varying at large radii, r > f20o/3 
(D99), when assuming the NFW profile for the cluster 
mass density and computed the anisotropy profile j3{r) 
from simulations. D99 followed DG97 and set to be 
0.5, finding that the resulting caustic method recovers 
the actual cluster mass within a factor of 2 at large radii, 
r ~ (0.3 — 6) r2oo, at least with the low concentration pa- 
rameters of his simulated clusters. Figure [8] shows that 
for A1689 there is a reasonably good agreement between 
the mass profile derived from the caustics and that de- 
rived from our earlier joint lensing and X-ray analysis, 
except (for the uncorrected caustic method) in the inner 
region, although the data points even in this region are 
within the 1-sigma range of uncertainty. The difference 
is well below a factor of 2 in the range ~ (0.1 — 1.5) h^^ 
Mpc. This suggests that the assumption of Fp = 0.5 
is reasonable for this cluster over a broad radial range. 
Note that the NFW profile fitted to this cluster gives a 
relatively high concentration parameter (see §[3]). 

At smaller radii, r < 0.1 h~^ Mpc, the two mass pro- 
files differ by ^ 1-a. This difference can arise if F^ is 
significantly different from 0.5. As a consistency check, 
we can use the velocity anisotropy profile and the mass 
profile of the cluster that we derived from the lensing 
and X-ray analysis to calculate the profile. Figure [5] 
shows that at small radii (r < 0.1 h"-'^ Mpc) Fp is indeed 
well below 0.5. In Figure [5] we compare the FJg profile 
derived by our model-independent method to that de- 
rived from assuming an NFW profile for the total mass 
density. The profile determined by D99 from simulations 
is somewhere in between our two measured profiles: the 
simulated F/^ increases with radius and then flattens, but 
only at large radii (~ j^-vir). Note that in the model- 
independent method, the point at the largest radius is 
lower than in the NFW case, since the density there is 
lower than derived using the NFW model due to the fact 
that the slope of the total density at large radii for A1689 
is slightly steeper than for the NFW profile (B05a; LOS). 
In Figure [S] we also plot the mass profile based on the 
caustics method when we use the more realistic Fp value 




10 

Radius [h^^ kpc] 

Fig. 8. — Profile of total enclosed mass, as derived from the 
caustics (red curve), caustics with corrected for r < 250 h~^ 
kpc (green dash-dotted curve), lensing and X-ray data (blue circles, 
LOS), and an NFW mass profile fit to the galaxy surface number 
density and projected velocity dispersion data sets (black dots). 
The uncertainties (error bars for the lensing and X-ray data based 
points, and red dashed curves for the caustic based profile) are l-cr. 



(instead of 0.5) at r < 250 h^^ kpc. The correction im- 
proves the agreement between the caustics-based mass 
profile and the profile derived from the joint lensing and 
X-ray analysis of L08. 




100 

Radius [h^^ l^pc] 

Fig. 9. — Reconstructed profile of the dimensionless factor Ffj 
given by eq. lO. We compare two profiles, where we use the total 
mass density derived by the model-independent method in LOS 
(blue points) or the best-fit NFW profile from LOS (red curve). The 
uncertainties (blue error bars for the model-independent method 
and dashed curves for the NFW model) are l-cr. Also shown for 
comparison is the F/j = 0.5 (dotted horizontal) line. 

The mass of the cluster derived by the caustics method 
is roughly equal to the virial mass, Afvir = (1.23±0.42) x 
10^^ h"-'^ M0, since the caustics end at about the virial 
radius (see §[6]). 

5.2. Mass profile from the Jeans equation 

Analysis of the lensing data suggests that A1689 has a 
high concentration parameter (B05a; LOS). We can in- 
dependently check the concentration parameter and the 
virial mass as derived solely from our two data sets that 
measure the galaxy surface number density and projected 
velocity dispersion. A similar attempt to do so was made 
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by Lokas et al. (2006) but with significantly less exten- 
sive projected velocity data (see §2.2|) . In addition to 
what has been done in § [31 we have also obtained the 
mass by implementing the following procedure: We use 
the galaxy number density and velocity anisotropy pro- 
files (eqs.[6]and[7]), along with an NFW profile for M{r), 
which adds free parameters Cvir (the concentration pa- 
rameter) and Mvir (the virial mass), and derive the pa- 
rameters of the three profiles by simultaneously fitting to 
the galaxy number density and velocity data, using the 
Jeans equation. This method results in relatively weak 
constraints, as is clear from the rather wide 1-cr contours 
shown in figure [10] (black contours); the deduced mass 
and concentration are Mvir = l-6lj g x 10^^ h~^, and 
Cvir > 5.8 (with no useful upper limit). 




1 1.5 2 2.5 3 3.5 4 

Mass IIO'^ h"' M ] 



Fig. 10. — Virial mass and concentration parameter, Cvir, derived 
from NFW fits to the cluster total mass profile. We show the results 
of an NFW fit to the galaxy dynamics using the Jeans equation 
(black circle for the best-fit values and black contour for the 1-cr 
uncertainty), the mass measured at the virial radius from caustics 
(red solid and dashed lines for the best-fit value and l-cr upper 
limit), and the joint 1-cr uncertainty taking both methods, using 
the Jeans equation and the caustics, under consideration (green 
contours). For comparison we show the result from the lensing 
and X-ray data (blue square with 1-cr error bars; LOS). 

The latter method of estimating the cluster mass, from 
the Jeans equation with assumed profiles, is less precise 
than the other methods, yielding a partial degeneracy 
between Cvir and Mvir as shown in figure [TUl We thus ex- 
pect this constraint to be weak. The observed galaxy 
number density essentially determines ?T.gai('') through 
eq. ([U, and the observed projected velocity dispersion 
determines a degenerate combination of and p at each 
radius through eq. This yields a degeneracy where 
for any assumed /3{r), the Jeans equation yields an M(r) 
that is consistent with the galaxy dynamical data. In the 
actual fitting though, this degeneracy is partially broken 
by the strict analytical forms assumed for the various in- 
put profiles. Figure [5] shows that the best-fit NFW mass 
profile from the Jeans equation is in good agreement with 
the mass profiles from the other methods. 

To get more useful constraints on Cvir and Afvir, we 
combine the two dynamical methods, i.e., the caustics 
and the Jeans equation. From the caustics method we 
take only the constraint on Mvir, since the mass pro- 
file at small radii is uncertain in this method due to the 
breakdown of the Ffj = 0.5 assumption (see discussion in 



the previous subsection). Figure [TOl shows that the com- 
bined constraints (green contours) are stronger and in 
good agreement with the values derived from the lensing 
and X-ray data by LOS (who assumed in this particu- 
lar analysis an NFW profile for the total mass density 
and a double beta model for the gas mass density pro- 
file). In particular, the combined dynamical methods 
yield 1-cr limits of Mvir = (1.3 ±0.4) x 10^^ h'^ M© and 
Cvir > 13.4. 

6. THE EDGE OF THE CLUSTER 

In this section we summarize various ways of defin- 
ing a limiting radius for A1689. One way to define the 
edge of the cluster is to use the observed galaxy number 
density profile. The fits that we have used in eqs. ([l} 
and ([6]) do not have a sharp cutoff in the number den- 
sity of galaxies. We thus define an edge as the radius 
where we can no longer detect cluster members above 
the contribution of the (non-cluster) background galaxy 
level. For A1689 this point is visible in figure [1] Specif- 
ically, the cluster radial edge was estimated to be where 
Sgai/Cbg > ACbg/Cbg and Sgai/Cbg < ACbg/Cbg in- 
cluding uncertainties in Sgai/Cbg, yielding a limiting ra- 
dius of 2.1^q'7 h~^ Mpc, where the 1-cr uncertainties ac- 
count also for the errors in the various fitting parameters 
in eq. 

Independently, the velocity caustic fits to the projected 
velocity dispersion data shown in figure [2] yield a very 
similar value for the limiting radius, 2.12±0.07 h~^ Mpc, 
where the error includes an estimate of the effect of Pois- 
son noise in the observed number of galaxies. We caution 
that in simulations the caustics often fiatten but do not 
reach zero at the virial radius; also, the shape of the 
caustics is somewhat dependent on the particular line of 
sight (D99). However, the caustics are generally more 
cleanly defined in data on real clusters than in N-body 
simulations (Rines et al. 2003). 

Both of these methods yield a cluster edge limiting 
radius of ^ 2 h~^ Mpc. A similar value was also inde- 
pendently derived from our lensing and X-ray analysis, 
which depends mostly on the projected DM distribution. 
We found the virial radius to be 2.14+";29 h'^ Mpc (L08). 
We conclude that all these different data sets agree rea- 
sonably well both in terms of the virial mass and mass 
profile and in the values they yield for the cluster's lim- 
iting radius. 

The various ways of defining a limiting radius for 
A1689 are consistent with each other. This suggests that 
there is no major infall of DM and galaxies. Significant 
ongoing infall would add to the projected profiles of mass 
and galaxy number and also affect the dynamical mea- 
surements, likely making the cluster edge less apparent. 

7. DISCUSSION 

In this paper we have continued our exploration of 
A1689 making use of many high quality datasets avail- 
able for this cluster. This work builds upon our earlier 
work on this cluster (L08), where we developed a com- 
prehensive joint analysis of high quality strong lensing 
(HST/ACS), weak lensing (Subaru), and X-ray (Chan- 
dra) measurements, from which we tested the consistency 
of X-ray and lensing data in a model independent way 
and derived an improved mass profile for A1689. In this 
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paper we have incorporated two other high quahty data 
sets, the galaxy surface number density measured from 
deep, wide-field imaging with Subaru/Suprime-Cam and 
a large spectroscopic study of the internal galaxy dy- 
namics measured using VLT/VIMOS. While the lousing 
and X-ray data gave us the information on the DM and 
gas content of the cluster, the two new data sets added 
here provide direct information on the galaxy distribu- 
tion and motions, leading to new determinations of the 
DM distribution of the cluster. 

The 3D galaxy number density profile derived from 
our combined analysis of the above datasets is more con- 
sistent with a cored profile, rather than a cuspy profile 
such as the NFW profile, which seems to fit well the DM 
density distribution. This is in agreement with Adami 
et al. (1998) who examined a sample of 62 clusters and 
found that the majority are better fit with a core than 
a cuspy profile, though for individual clusters the pref- 
erence for a cored profile is rarely significant at the 90% 
confidence level. The galaxy distribution also resembles 
a King (1962) profile and falls off as 7.-3 i8±o.42^ exceed- 
ing the slope of suggested by Bahcall & Lubin 
(1994) in order to explain the "/3-discrepancy" . This 
asymptotic behavior of the cluster galaxy profile at large 
r is in fact very similar to our total matter profile (dom- 
inated by DM), which we have shown can be well fit by 
an NFW profile where the asymptotic behavior is . 

A principal finding of our work is the first direct deter- 
mination of the velocity anisotropy profile for a galaxy 
cluster. This followed from an application of the Jeans 
equation, using as input the observed projected velocity 
dispersion profile, the observed projected galaxy distri- 
bution, and our independently determined mass profile, 
allowing us to solve for the 3D velocity anisotropy as 
a function of radius. The resulting anisotropy profile 
is well fit by the expression in equation ([7]) proposed 
by Carlberg et al. (1997) on the basis of N-body sim- 
ulations. The simulations covered a wide range of cos- 
mologies and showed that the radial dependence of the 
velocity anisotropy (3 has a nearly universal form (Cole 
& Lacey 1996; Carlberg et al. 1997), with a characteris- 
tic radial dependence. This dependence is also observed 
in A1689 - mainly radial motion is deduced at large ra- 
dius, tending towards isotropic (or possibly tangential) 
motion within the central region r < 500 h~^ kpc. This 
presumably is a manifestation of the overall formation 
and growth of clusters, with initial collapse and virial- 
ization of the central region, and continued growth of 
the cluster mass through accretion, a two-stage process 
for which there seems to be some evidence also from the 
DM entropy distribution (e.g., Lapi & Cavaliere 2008). 

Additionally, we used new extensive measurements of 
galaxy positions and velocities to determine the cluster 
mass profile. This was done in two independent ways; 
first, using the velocity data alone, we identified clearly 
apparent velocity caustics using the method of D99. The 
derived amplitude of the velocity caustics was interpreted 
as the local escape velocity, from which the mass profile 
was determined using eq. Secondly, we followed the 
traditional approach of using the Jeans equation (eq. [2]) , 
incorporating both the galaxy surface number density 
and the projected velocity data, and adopting the above 
velocity anisotropy profile. These two different mass es- 
timates are in good agreement with the profile derived 



from our earlier lensing and X-ray analysis (L08), as 
shown in figure [51 

In estimating the caustic-based mass profile we 
adopted the previously suggested value Fp = 0.5. We 
were able to separately check this assumption using the 
velocity anisotropy profile obtained as described above. 
Our previously determined mass profile made it possible 
to deduce /?; the comparison of the resulting mass pro- 
file with that from our previous lensing/X-ray analysis is 
then essentially a consistency check on the general valid- 
ity of the caustic method. We found that at large radii, 
r > 100 h-i kpc, the factor Ffs varies slowly with radius 
and stays within ~ 50% of the value of 0.5 (see figured]), 
so that the relatively simple caustic method yields the 
mass profile accurately except at the center. Note that 
D99 derived the mean value of Fp = 0.5 from simulated 
CDM halos, which typically have concentration parame- 
ter < 7, well below that of A1689 (cvir = 12.2;?-^; L08). 
Thus, our results suggest that the caustic-based mass 
estimation is applicable also for high concentration clus- 
ters. This shows that dynamical analysis may be im- 
proved upon by combining the traditional method based 
on the Jeans equation with the additional insight gained 
from the caustics (in relaxed clusters). 

The virial mass derived from the caustic method, 
Mvir = (1.23 ±0.42) X 10^5 h-i Mq, and from the Jeans 
equation, A/vir = l-Gto,l ^ 10^^ ^0- The two com- 
bined dynamical methods gave Afvir — (1.3 ± 0.4) x 10^^ 
h^^ Mq, as compared with the value obtained in our joint 
lensing/X-ray analysis (L08), Afvh- = (1.58 ±0.15) x lO^^ 
h~^ Mq. These mass estimates are consistent, and in 
agreement also with the result of Umetsu & Broad- 
hurst (2008), who combined strong lensing, weak lens- 
ing distortion and magnification data in a 2D anal- 
ysis (without assuming axial symmetry), and derived 
M^i, = l.5toi X 10^^ Mq (where this 1 - cr error 
includes both statistical and systematic uncertainties). 

A novel aspect of our work is an estimation of the lim- 
iting radius of A1689 in several different ways. All our 
data sets independently indicate that this cluster has a 
relatively well defined edge radius of ~ 2.1 h"-*^ Mpc. We 
see this in terms of the galaxy distribution which is flat 
beyond this radius, indicating that the galaxy number 
density reached the background level, from the velocity 
caustics which attain (practically) zero velocity at this 
radius, and also from the virial radius derived from the 
best-fit NFW profile. This interesting agreement should 
be checked in other massive clusters. In principle, we 
may expect a clear signature of the virial radius in the 
absence of significant ongoing infall of galaxies and mass 
onto the cluster. 

Massive clusters {M ^ 10^^ h~^ Mq) are theoretically 
expected to have relatively low mass concentrations, with 
Cvir(z = 0) ~ 5, as determined from cosmological simula- 
tions of the standard ACDM (Neto et al. 2007; Hennawi 
et al. 2007; Duffy et al. 2008). However, significantly 
higher values, Cvir ~ 10—15, were deduced from de- 
tailed weak and strong lensing measurements of several 
well-studied massive clusters (Kneib et al. 2003; Gavazzi 
et al. 2003; B05a; Kling et al. 2005; Limousin et al. 
2007; L08; Umetsu & Broadhurst 2008; Broadhurst et al. 
2008). This discrepancy can also be expressed in terms of 
the Einstein radius, which provides a model-independent 
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measure of the central mass within clusters, and can be 
compared with the total (virial) mass to get a measure 
of the degree of concentration. The observed Einstein 
radii of many massive clusters lie in the range 20" — 50", 
around twice the expected range after allowance for lens- 
ing and projection biases (Broadhurst feBarkana 2008). 
In addition, Sadeh & Rephaeli (2008) showed that even 
considering the full probability distribution function of 
halo formation times does not fully remove the discrep- 
ancy. 

While we have not measured the concentration param- 
eter of A1689 as accurately as from the Icnsing data, we 
have confirmed its high value. The two combined dy- 
namical methods yield cvir > 13.4 (at 1-a confidence), 
consistent with c^i^ = 12.2^q g derived in our joint lens- 
ing and X-ray analysis (L08), and with Cvir = 12.7 ± 3 
from Umetsu & Broadhurst (2008). Clearly, this com- 
parison should be examined for other well-studied clus- 
ters, but at least in A1689 this independent dynamical 
measurement, which is less susceptible to projection bias 
than lensing, is in good agreement with the lensing re- 
sults, and can be seen (at least partly) as evidence for 
the validity of (the assumption of) spherical symmetry. 

The effect of triaxiality may significantly lower the ob- 
servationally deduced concentration parameter (Oguri et 
al. 2005; Gavazzi 2005). An independent analysis of 
A1689 by Corless, King, & Clowe (2008) does not find 
this to be significant, Cmr = 12.2 ± 6.7, in good agree- 
ment with our earlier analysis, where the sizable uncer- 
tainty reflects the unknown triaxiality and the relatively 
shallow imaging used. An important step in empirically 
examining the effect of triaxiality and geometric projec- 
tion bias is to see whether there is any systematic dif- 
ference between the distribution of mass profiles derived 
for clusters selected by different methods. In a recent 
analysis of a cluster sample selected rising weak lensing, 
Hamana et al. (2008) who used the sample of Miyazaki 
et al. (2002) - did not find evidence for selection bias in 
mass estimates between clusters selected by weak lens- 
ing as compared with other selection methods, such as 
optically or by X-ray emission. Also, Duffy et al. (2008) 
found that the derived concentrations for X-ray selected 
samples taken from the literature lie well above the ex- 
pected values based on cosmological N-body simulations 
with the most recent cosmological parameters. 

At least some of the high deduced values of Cmr are ex- 
pected in early dark energy (EDE) models. In these mod- 
els there is a non- negligible dark energy component even 
at high redshift, in contrast with the standard model; 
consequently, clusters form earlier and are more abun- 
dant than predicted in ACDM (e.g, Sadeh & Rephaeli 
2008). In recent work using N-body simulations it was 
found that halos in an EDE model show higher concen- 
trations at a given halo mass (Grossi & Springel 2008). 
Interestingly, in contrast with expectations from semi- 
analytic approaches, recent simulations show that EDE 
does not significantly affect the statistics of non-linear 
structures, at least at low redshifts (Francis, Lewis, & 
Linder 2008; Grossi & Springel 2008). In addition, the 
measured X-ray background cannot rule out EDE mod- 
els if their normalization is defined by the halo density 
at z = (Lemzc, Sadeh, & Rephaeli 2009). We intend 
to explore further the issue of values of Cyir by extending 



our dynamical work to additional well-studied clusters 
for which we have obtained precisely measured lensing- 
based mass profiles. 
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